NAVAL 

POSTGRADUATE 

SCHOOL 

MONTEREY,  CALIFORNIA 


THESIS 


LAGUERRE-GAUSSIAN  MODES 
IN  THE  FREE  ELECTRON  LASER 

by 

Anastasios  Kampouridis 
June  2007 

Thesis  Advisor:  William  B.  Colson 

Co-Advisor:  Robert  L.  Armstead 


Approved  for  public  release;  distribution  is  unlimited 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


REPORT  DOCUMENTATION  PAGE 


FomiApprovedO/WBWa_0704^0f88_ 
Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing 
instruction,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection 
of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden,  to  Washington  headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215 
Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction 
Project  (0704-0188)  Washington  DC  20503. 

1 .  AGENCY  USE  ONLY  (Leave  blank)  I  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

June  2007  Master’s  Thesis 

4.  TITLE  AND  SUBTITLE  Laguerre-Gaussian  Modes  in  the  Free  Electron  5.  FUNDING  NUMBERS 
Laser _ 

6.  AUTHOR(S)  Anastasios  Kampouridis _ 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Naval  Postgraduate  School  REPORT  NUMBER 

Monterey,  CA  93943-5000 

9.  SPONSORING  /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSORING/MONITORING 
N/A  AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the 
official  policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. _ 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE  = 

Approved  for  public  release;  distribution  is  unlimited _ 

13.  ABSTRACT  (maximum  200  words) 

In  a  free  electron  laser  (FEL)  system,  knowing  the  optical  beam  characteristics  is  of  great  importance.  A  beam 
may  be  comprised  of  higher-order  modes  due  to  the  interaction  with  the  electron  beam,  or  from  non-ideal  operational 
conditions  such  as  mirror  distortions  and  misalignments,  or  from  imperfect  injection  of  the  electron  beam. 

In  this  thesis,  the  basic  FEL  theory  is  initially  reviewed.  The  parabolic  wave  equation  is  then  solved  for  the 
“fundamental”  Gaussian  mode  and  for  higher-order  modes.  Working  in  rectangular  coordinates,  a  complete  and 
orthogonal  set  of  solutions  involving  Hermite  polynomials  is  found.  When  the  wave  equation  is  solved  in  cylindrical 
coordinates,  we  arrive  at  a  set  of  solutions  that  contain  Laguerre  polynomials.  The  so-called  Laguerre-Gaussian  modes 
are  analyzed.  The  evolution  of  these  laser  modes  is  also  explored,  yielding  quite  unexpected  results  due  to  their  phase 
structure  and  orbital  angular  momentum  of  light.  Lastly,  we  study  a  common  case  where  higher-order  optical  modes 
appear,  in  order  to  quantify  the  tolerances  of  an  FEL. 


14.  SUBJECT  TERMS  Free  Electron  Laser,  High  Order  Modes,  Laguerre-Gaussian  Optical  15.  NUMBER  OF 
Modes.  PAGES 

87 

16.  PRICE  CODE 

17.  SECURITY  18.  SECURITY  19.  SECURITY  20.  LIMITATION  OF 

CLASSIFICATION  OF  CLASSIFICATION  OF  THIS  CLASSIFICATION  OF  ABSTRACT 

REPORT  PAGE  ABSTRACT 

_ Unclassified _ Unclassified _ Unclassified _ UL _ 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  239-18 


I 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


Approved  for  public  release;  distribution  is  unlimited 


Author: 


Approved  by: 


LAGUERRE-GAUSSIAN  MODES 
IN  THE  FREE  ELECTRON  LASER 

Anastasios  Kampouridis 
Major,  Hellenic  Air  Force 
B.S.,  Hellenic  Air  Force  Academy,  1991 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  PHYSICS 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
June  2007 


Anastasios  Kampouridis 


William  B.  Colson 
Thesis  Advisor 


Robert  L.  Armstead 
Co-Advisor 


James  H.  Luscombe 
Chairman,  Department  of  Physics 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


IV 


ABSTRACT 


In  a  free  electron  laser  (FEL)  system,  knowing  the  optical  beam 
characteristics  is  of  great  importance.  A  beam  may  be  comprised  of  higher-order 
modes  due  to  the  interaction  with  the  electron  beam,  or  from  non-ideal 
operational  conditions  such  as  mirror  distortions  and  misalignments,  or  from 
imperfect  injection  of  the  electron  beam. 

In  this  thesis,  the  basic  FEL  theory  is  initially  reviewed.  The  parabolic 
wave  equation  is  then  solved  for  the  “fundamental”  Gaussian  mode  and  for 
higher-order  modes.  Working  in  rectangular  coordinates,  a  complete  and 
orthogonal  set  of  solutions  involving  Hermite  polynomials  is  found.  When  the 
wave  equation  is  solved  in  cylindrical  coordinates,  we  arrive  at  a  set  of  solutions 
that  contain  Laguerre  polynomials.  The  so-called  Laguerre-Gaussian  modes  are 
analyzed.  The  evolution  of  these  laser  modes  is  also  explored,  yielding  quite 
unexpected  results  due  to  their  phase  structure  and  orbital  angular  momentum  of 
light.  Lastly,  we  study  a  common  case  where  higher-order  optical  modes  appear, 
in  order  to  quantify  the  tolerances  of  an  FEL. 
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I.  INTRODUCTION 


The  Free  Electron  Laser  (FEL)  is  a  source  of  powerful,  coherent  radiation 
with  unique  attributes.  Continuous  tunability,  high  power,  good  efficiency  and 
reliability  make  it  suitable  for  a  wide  span  of  applications.  The  option  of  high 
power  derives  in  part  from  the  absence  of  a  medium  vulnerable  to  damage. 

Powerful  coherent  radiation  was  first  produced  when  the  microwave  tube, 
which  uses  a  beam  of  free  non-relativistic  electrons  inside  a  closed  microwave 
cavity,  was  invented  in  the  1930’s.  Following  this  invention,  the  introduction  of  the 
open  resonator  enabled  many  steps  towards  the  development  of  more  powerful 
and  efficient  radiation  sources. 

Based  on  its  special  attributes,  the  FEL  is  suitable  for  many  diverse 
applications.  It  can  be  used  in  biology  to  study  cell  and  molecular  structure,  and  it 
can  be  used  in  the  military,  due  to  its  high  power  output  and  its  capability  to 
operate  in  a  broad  range  of  wavelengths.  A  critical  requirement  for  an  FEL  to 
have  the  desired  results  is  good  optical  beam  quality.  We  need  to  know  the  beam 
shape  that  can  be  produced  under  various  operational  conditions,  so  as  to 
evaluate  its  performance. 

In  this  thesis,  we  analyze  various  high-order  modes  from  an  FEL, 
emphasizing  cylindrical  coordinates,  which  are  motivated  by  the  symmetry  of  the 
system.  We  present  a  brief  review  of  a  typical  FEL  system  in  the  first  chapter,  as 
well  as  the  relevant  operating  configurations.  The  next  chapter  introduces  basic 
FEL  theory,  where  we  derive  the  pendulum  equation,  which  describes  the 
electron  motion,  and  the  parabolic  wave  equation,  which  describes  the  laser  light 
evolution,  from  the  classical  electromagnetic  theory. 

In  the  chapters  that  follow,  we  focus  our  attention  on  the  optical  beam. 
First,  the  basic  optical  theory  is  reviewed  and  the  fundamental-mode  solution  to 
the  wave  equation  is  derived  and  analyzed.  Although  this  solution  is  a  good 
approximation  for  the  output  of  many  FELs,  the  actual  operating  parameters  of  a 
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system  can  generate  laser  light  with  higher-order  modes.  To  analyze  this  effect, 
we  initially  review  work  with  rectangular  coordinates,  and  find  the  Hermite- 
Gaussian  modes.  Following  this,  new  work  is  presented  with  cylindrical 
coordinates,  which  are  more  relevant  to  an  FEL  system,  and  we  arrive  at  the 
Laguerre-Gaussian  solutions.  Analyzing  these  high-order  modes,  we  find  that  the 
orbital  angular  momentum  of  light  plays  important  role  in  the  evolution  of  such  a 
beam. 

In  the  last  chapter,  we  present  the  results  of  a  new  set  of  simulations 
designed  to  investigate  the  tolerance  of  the  Jefferson  Lab  FEL  (Figure  1)  to 
electron  beam  tilt.  From  these  results,  we  see  that  higher-order  optical  modes  do 
indeed  appear  as  a  consequence  of  the  non-ideal  operating  conditions. 


Electron  Gun  High  voltage 

power  supply 


Figure  1 .  Schematic  of  Jefferson  Lab  free  electron  laser  After  [1], 
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II.  FREE  ELECTRON  LASER  SYSTEM 


In  the  Free  Electron  Laser,  a  beam  of  relativistic  electrons,  produced  by 
an  injector  and  accelerated  by  a  linear  electron  accelerator,  passes  through  a 
transverse,  periodic  magnetic  field,  called  the  undulator,  oscillates  and  therefore 
radiates.  Radiation  can  be  captured  in  a  cavity  and  used  to  induce  new  electrons 
to  emit  even  more  radiation.  Alternatively,  radiation  from  an  external  source  can 
be  amplified  over  a  single  pass,  and  even  without  an  external  source  it  is 
possible  to  produce  radiation  through  the  process  of  spontaneous  emission. 
Therefore,  it  is  informative  to  look  at  the  possible  basic  setups,  which  can  be 
used  in  order  for  a  free  electron  laser  system  to  function. 

A.  FEL  CONFIGURATIONS 

The  configuration  used  in  the  first  FEL  by  John  Madey  in  1977  is  the 
oscillator,  employing  an  optical  cavity  in  a  very  similar  way  as  in  the  conventional 
chemical  or  solid  state  lasers.  Radiation  is  stored  and  amplified  in  the  optical 
cavity,  bracketed  by  one  fully  and  one  partially  reflective  mirror,  the  latter  being 
used  to  outcouple  the  laser  light.  Coherence  is  established  after  many  passes  of 
the  light.  Figure  2  is  a  typical  FEL  oscillator,  where  the  emphasis  is  placed  on  the 
resonator. 


mirror 


electron  source 
(gun) 


Figure  2.  Oscillator  Free  Electron  Laser  From  [2], 
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In  the  amplifier  configuration  there  is  no  optical  resonator  and  it  can  be 
considered  the  simplest  type  of  FEL.  Light  from  a  seed  laser  is  amplified  by  the 
interaction  with  the  electron  beam,  passing  only  once  through  the  periodic 
magnetic  field.  Coherence  in  this  case  is  set  up  by  the  seed  laser. 

A  special  type  of  amplifier  is  the  SASE  FEL  (Self-Amplified  Spontaneous 
Emission)  where  there  is  no  external  source.  Instead,  the  spontaneous  radiation 
in  the  first  part  of  the  magnetic  field  is  used  as  the  seed  laser.  At  present,  It 
appears  to  be  the  most  promising  concept  in  the  race  towards  achieving  much 
shorter  wavelengths  (UV, X-rays),  since  there  are  no  seed  lasers  in  this  regime. 
Figure  3  shows  the  layout  of  a  SASE  FEL,  where  there  are  neither  resonator 
mirrors  nor  seed  laser,  but  where  a  significantly  longer  undulator  is  required. 


Figure  3.  SASE  Free  Electron  Laser  From  [2], 

B.  FEL  COMPONENTS 

The  major  components  of  an  FEL  system  are  the  electron  beam  source, 
the  accelerator,  the  undulator,  the  resonator  (for  the  oscillator)  and  the  beam 
dump.  There  are  many  other  elements  that  make  an  FEL  system  operational,  like 
bending  and  control  magnets  to  steer  the  electron  beam  in  the  desired  path  and 
several  detecting  devices  to  continuously  monitor  all  the  critical  points,  in  order  to 
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keep  track  of  the  operational  parameters  and  to  enhance  the  performance.  Last 
but  not  least,  an  optical  system  must  be  utilized  to  transport  the  laser  light  in  the 
desired  area. 

1.  Electron  Beam  Injector 

The  electron  beam  injector  is  the  original  source  of  electrons  in  the  FEL, 
and  it  plays  a  crucial  roll  in  the  performance  of  the  system.  It  is  usually  a 
photoinjector,  where  a  drive  laser  excites  the  electrons  in  a  cathode,  which  is 
positioned  inside  a  Radio  Frequency  (RF)  cavity,  and  a  high  voltage  is  applied  so 
as  to  accelerate  the  released  electrons,  allowing  them  to  gain  a  modest  amount 
of  energy  before  entering  the  main  accelerator.  Instead  of  taking  advantage  of 
the  photoelectric  effect,  much  the  same  release  of  electrons  can  be  obtained  by 
heating  the  cathode.  The  electron  beam  that  is  produced  from  an  injector,  like  the 
one  shown  in  Figure  4,  usually  comes  out  in  pulses  at  regular  intervals. 


Figure  4.  Electron  Injector  After  [3], 

2.  Accelerator 

The  electrons  emerge  from  the  injector  with  energy  of  about  5-10  MeV.  In 

order  to  achieve  infrared  radiation  from  the  FEL,  we  need  a  relativistic  electron 

beam  of  about  100  MeV.  To  attain  this  energy  level,  the  electron  beam  is  passed 

through  a  Radio  Frequency  Linear  Accelerator  (RF  LINAC).  Superconducting 
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technology  may  be  used  as  in  the  JLab  FEL  LINAC,  shown  in  Figure  5.  The 
accelerator  uses  an  alternating  electric  field,  which  interacts  with  each  electron 
bunch  that  passes  through  the  RF  cavity.  For  this  reason,  the  electric  field  has  to 
be  fully  synchronized  with  the  electron  bunches,  so  that  it  continuously 
accelerates  the  electrons.  Also,  the  accelerator  must  have  several  stages  of  RF 
modules  in  order  for  the  electrons  to  reach  the  required  energy. 

The  RF  field  inside  the  linear  accelerator  is  generated  by  a  klystron.  To 
reduce  the  input  power  required,  a  special  technique  has  been  introduced.  The 
electron  beam  that  has  traveled  through  the  whole  system  can  be  driven  back 
through  the  accelerator,  so  as  to  give  most  of  its  remaining  energy,  back  to  the 
RF  field.  To  achieve  this,  we  have  to  carefully  adjust  the  timing  of  the  beam,  so 
as  to  be  out  of  phase  with  the  electric  field  of  the  RF  cavities.  This  energy 
recovery  procedure  through  the  recirculation  of  the  electron  beam  increases  the 
overall  efficiency  of  the  system  and  additionally  reduces  the  beam  dump  size. 


(a) 


(c) 

Figure  5.  JLab  Superconducting  LINAC  (a),  (b)  A  five-cell  cavity,  photo  and 
sketch,  (c)  Part  of  a  cryomodule  which  contains  six  five-cell  cavities 
From  [4], 
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3. 


Undulator 


The  undulator  is  the  component  of  the  FEL  system  that  forces  the  electron 
beam  to  wiggle  and  radiate.  It  is  composed  of  a  group  of  magnets,  as  can  be 
seen  in  Figures  2  and  3,  arranged  so  that  they  form  an  alternating  magnetic  field 
along  the  axis  of  the  undulator.  The  electrons  leaving  the  accelerator  with 
relativistic  speeds,  enter  the  periodic  magnetic  field  of  the  undulator,  where  the 
Lorentz  force  drives  them  to  oscillate  in  a  transverse  direction.  The  radiation  that 
is  emitted  interacts  with  the  light  already  present  in  the  undulator,  through  the 
process  of  the  stimulated  emission.  This  interaction  results  in  the  amplification  of 
the  stored  light. 

4.  Resonator 

In  the  oscillator  configuration,  as  the  one  shown  in  Figure  2,  the  FEL 
system  makes  use  of  an  optical  cavity  with  the  undulator  enclosed  inside.  This 
way,  light  is  stored  and  amplified  through  many  passes.  Two  mirrors  bound  the 
cavity:  one  of  them  totally  reflects  the  light,  while  the  other  allows  a  percentage  to 
come  out.  The  light  that  emerges  through  this  out-coupling  mirror  is  the  desired, 
usable  laser  light  of  the  system.  It  is  obvious  that  for  maximum  coupling 
efficiency  between  the  electrons  and  the  light  within  the  resonator,  the  optical 
beam  should  be  synchronized  with  the  electron  bunch. 

5.  Beam  Dump 

The  beam  dump  is  simply  a  notched  piece  of  metal  that  absorbs  the 
energy  of  the  electrons  at  the  end  of  their  path,  either  just  after  the  undulator  or  in 
the  energy  recovery  configuration,  after  the  linear  accelerator.  Even  though 
electrons  give  energy  to  the  optical  field  in  order  for  the  FEL  to  emit  light,  and 
probably  back  to  the  RF  field,  they  still  have  energy  to  be  removed.  Therefore, 
the  beam  dump  can  heat  up  at  megawatt  rates  and  needs  to  be  water-cooled. 
Furthermore,  it  is  shielded  so  that  no  radiation  can  exit.  A  schematic  diagram  of  a 
beam  dump  is  shown  in  Figure  6. 
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Schematic  diagram  of  the  electron  beam  dump  After  [5], 


8 


III.  BASIC  FREE  ELECTRON  LASER  THEORY 


The  basic  idea  behind  the  operation  of  a  free  electron  laser  is  the 
extraction  of  energy  from  relativistic  electrons,  forcing  them  to  radiate  coherently. 
This  energy  conversion  takes  place  inside  the  undulator,  where  the  electron 
beam  makes  a  wiggling,  transverse  motion,  emitting  spontaneous  radiation  in  the 
forward  direction.  We  can  couple  this  motion  to  the  electric  field  of  the  radiation 
to  produce  stimulated  emission,  thus  amplifying  the  optical  mode.  We  will  now 
look  at  the  mathematics  that  describe  the  attributes  of  a  free  electron  laser  and 
explain  its  operation. 

A.  FEL  PENDULUM  EQUATION 

In  order  to  understand  the  physics  behind  the  electron  motion  and  the 
subsequent  radiation,  we  will  first  study  their  behavior  inside  the  magnetic  field  of 
the  undulator.  For  that  reason,  we  will  consider  relativistic  electrons  entering  a 
helical  undulator  field  along  the  z-axis  in  the  presence  of  optical  field.  As  the 
electrons  travel  along  the  undulator,  forces  from  both  the  static  magnetic  field  of 
the  undulator,  and  the  oscillating  electric  and  magnetic  fields  of  the  optical  mode 
affect  them  [6], 

The  magnetic  field  of  the  undulator  can  be  expressed  in  rectangular 
coordinates  as: 

By  =  B  (cos  (kaz  ),  sin  ( k0z ),  0  ),  ( 1 1 1 . 1 ) 

where  B  is  the  magnetic  field  strength,  k0  =2 k/Ao  is  the  undulator  wavenumber 
and  / L0  is  the  undulator  wavelength,  being  the  distance  of  one  complete  cycle  of 
the  magnetic  field  in  the  undulator. 

The  corresponding  electric  and  magnetic  field  of  the  optical  wave  are 
given  by: 
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E  =  E(cos  ^,-sin^-,0), 

(MI-2) 

B  =  E  (sin  y/,  cos  y/,  0  ) , 

(MI-3) 

with  V  -  kz  cot  +  cp ;  g  the  field  amplitude  of  the  optical  wave  in  cgs  units, 
k=oo!  c  the  optical  wavenumber,  co  the  frequency  and  <p  the  optical  phase. 

Electrons  entering  these  fields  encounter  a  Lorentz  force  that,  in  the  case 
of  relativistic  particles,  is  given  by  the  equations 


d{rP) 

dt 


e 

me 


(E  +  fixB), 


(HI-4) 


dy 

dt 


e 

me 


P-E 


(HI -5) 


where  y=  1 / ^l~p2  is  the  Lorentz  factor,  p  =  v/c  is  the  dimensionless 

electron  velocity,  e  =  \e\  is  the  electron  charge  magnitude  and  m  is  the  electron 

rest  mass.  Equation  (III. 5)  is  also  known  as  the  energy  equation  because  it 
describes  how  electrons  exchange  energy  with  the  optical  field,  and  the  electron 
energy  is  s  =  ymc 2 . 

We  can  now  substitute  the  fields  from  (III.  1  )-(lll.3)  into  (III .4)  to  obtain  the 
spatial  components  of  the  force  on  the  electron 


d(/P±) 

dt 


— —  [£(l-/?J(cos^,-sin^,0)  +  /LB( 
me 


sin£oz,cos&oz,0)] ,  (III. 6) 


dt 


e 

me 


E  (/?  cos  (//  -  Py  sin  i//j+  B  (j3x  sin  koz  -  Py  cos  k0z)  , 


(I  H-7) 
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where  Px,(3y,f3,  are  the  three  components  of  the  velocity  and  fi±  =  J3xi+  J3y  j  is 
the  transverse  velocity.  Since  we  are  dealing  with  relativistic  electrons,  we  can 
take  /?.  «  1 .  Then,  all  terms  of  E( 1- y 3 Q  are  small  compared  to  and  can  be 

dropped.  Integrating  (III. 6)  and  assuming  perfect  injection  into  helical  orbits,  so 
as  to  take  the  constants  of  integration  to  be  zero,  we  obtain  the  velocity  of  the 
transverse  motion 


Inymc2 


(cosA:oz,sinA:oz,0) . 


(II 1-8) 


We  can  write  (III. 8)  in  a  more  compact  form  if  we  introduce  the  dimensionless 
undulator  parameter 


K= 


eB . 1 


2  it  me 


(111.9) 


with  Bnm  the  rms  value  of  the  undulator  magnetic  field.  Then 


K_ 

Y 


(cosA:oz,sinA:oz,0) . 


(111.10) 


In  (111.10)  note  that  for  the  helical  undulator,  Bnns  =B  and  for  a  linearly  polarized 
undulator  Brms  =b/V2  . 

At  this  point,  we  are  able  to  look  at  the  basic  transverse  motion  of  the 
electrons,  by  making  the  assumption  that  they  travel  roughly  with  the  speed  of 
light  along  the  undulator.  Then  p:c~c  ,  and  taking  the  electron  velocity  to  be 

constant  we  can  write  z(r)« /3zct .  Using  a>o=  p:k0c~k0c  as  the  electron 
oscillation  frequency,  (111.10)  becomes 


—  K 

A~  (cos(<yor),sin(»or),0), 
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(111.11) 


and  integrating  we  get 


y  In 


(-sin(a>ot),cos(a>ot),0) , 


(111.12) 


with  xL(t)  =  x(t)i+y(t)j .  This  equation  approximates  the  transverse  motion  of 

the  electrons  inside  the  undulator.  From  this  expression  and  for  typical  systems, 
we  find  that  electrons  move  in  the  transverse  directions  in  a  distance  of  about 
25  jum ,  the  oscillation  frequency  is  in  the  microwave  regime,  but  due  to  the 

relativistic  Doppler  effect,  the  radiation  frequency  is  multiplied  by  a  factor  of  y2 , 
resulting  in  lasing  at  much  higher  infrared  frequencies. 

Next,  we  want  to  study  more  carefully  the  microscopic  motion  of  the 
electrons,  without  making  the  approximations  of  the  previous  paragraph.  For  this 
reason,  we  will  use  the  transverse  motion  we  found  in  equation  (111.10),  and  plug 
it  into  the  energy  equation  (III. 5)  to  get 

y  =  ^T=-— Efecos^-/?  sinyA  (111.13) 

dt  me  v  7 

Introducing  the  electron  phase  Q  =(k  +  ko^z  -  cot  and  recalling  that 
¥  =  kz  -  cot  +  rp ,  (|  1 1 . 1 3)  js  rewritten  as 

y=  (^^-cos(C+ <p\  (111.14) 

ymc 

The  last  equation  describes  energy  exchange,  because  it  tells  us  whether 
electrons  lose  or  gain  energy,  based  on  the  sign  of  the  cos(^+^)  term.  For 

instance,  when  cos(£+(p)<0  then  y<  0,  thus  electron  energy  s=ymc 1 

decreases,  which  is  the  desired  condition  for  the  FEL  to  operate,  so  that 
electrons  give  energy  to  the  optical  field. 
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Electron  phase  ^  is  a  way  to  relate  electron  position  to  the  co-propagating 
optical  field;  it  corresponds  to  the  position  of  the  electron  in  a  section  of  the 
electron  beam  one  optical  wavelength  long.  Therefore,  it  is  reasonable  to 
introduce  the  electron  phase  velocity,  defined  in  the  following  way  so  as  to  be 
dimensionless 


v  =  {  ±  =  L[(k  +  k')p._-k]  =  i ,  (111.15) 

where  L  =  NAa  is  the  length  of  the  undulator,  N  is  the  number  of  undulator 

periods  and  the  symbol  (..)  indicates  the  derivative  with  respect  to  the 
dimensionless  time  r,  defined  by  r  =  cr/Z  .  This  normalizes  time  to  the  length  of 
the  undulator,  so  that  r  =  0  -»  1  from  the  beginning  to  the  end  of  the  undulator. 

The  dimensionless  electron  phase  velocity  v  is  also  known  as  the 
“resonance  parameter”.  For  an  electron  that  satisfies  the  resonance  condition 
v=0,  the  electrons  and  optical  field  have  optimum  energy  exchange.  In  this 
case,  c/?  =  kc/(k  +  k0 ).  Assuming  again  relativistic  electrons,  we  have  /?z  «  1  and 
k^>kQ.  Therefore,  A  =  Aa  and  laser  frequency  ry  »  coa,  where  co0  «  2ncjAo  is  the 
electron  oscillation  frequency.  Then,  the  laser  wavelength  is 

A 


From  (111.10),  we  find  that  j3l=K2//2  and  so 


(111.17) 


For  relativistic  electrons  with  /  » 1 ,  equation  (111.17)  can  be  approximated  by 
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(111.18) 


A«i 


1 +  K2 


Combining  equations  (111.16)  and  (111.18)  we  find  the  following  expression  of  the 
resonance  condition  for  the  laser  wavelength 


X  = 


4(i+^2) 

If 


(111.19) 


This  resonance  relationship  makes  obvious  that  an  FEL  can  be  continuously 
tuned  by  changing  physical  parameters,  like  the  electron  energy  s  =  ymc 2  or  the 
design  of  the  undulator.  Since  the  undulator  parameter  is  typically  K  « 1 ,  it  is 
clear  that  optical  wavelength  is  much  smaller  than  that  of  the  undulator  by  a 
factor  of  approximately  \/y 2 . 

We  are  able  now  to  continue  our  analysis,  so  as  to  arrive  to  a  compact 
expression  for  the  electron  microscopic  motion,  taking  first  the  time  derivatives  of 
the  electron  phase  velocity  V,  defined  in  (111.15)  and  that  of  the  longitudinal 
velocity  /?_  from  (111.18) 


v  =  L(k  +  k0)pz , 


A  = 


1 +  K2 


~7  ■ 


(111.20) 


Then,  we  combine  these  results  to  find  v  =  4 xNy/y  and  finally  substituting  y 
from  (111.14)  we  get 


°  L  .  °°  AnNeKLE  , 

v  =  -v  =  £  = - ^ — —cos{<Z  +  (p). 

c  y  me 


(111.21) 


It  is  now  reasonable  to  collect  terms  and  define  |a|  =  47rNeKLE/y2mc2  as  the 

dimensionless  laser  field  amplitude,  so  that  (111.21 )  becomes 
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(111.22) 


O  00 

v=  |a|cos (£+<p), 

which  is  the  FEL  “pendulum  equation”  that  describes  the  microscopic  motion  of 
each  electron  in  the  reference  frame  of  one  optical  wavelength.  In  this  reference 
frame,  the  motion  of  the  electrons  can  be  taken  to  be  similar  to  a  classical 
pendulum.  The  dimensionless  laser  field  amplitude  |a|  is  a  measure  of  the 
strength  of  the  optical  field  and  expresses  the  rate  of  the  electron  bunching. 
When  |  a| »  k  the  laser  has  strong  fields  and  it  is  near  saturation,  while  for  |a|  <  n 
the  optical  field  is  weak  and  there  is  no  overbunching  of  electrons. 

Equation  (111.22)  implies  that  electrons  inside  the  magnetic  field  of  the 
undulator  move  in  an  equivalent  way  to  the  classical  pendulum.  Phase  space 
plots  are  commonly  utilized  to  visualize  the  motion  of  a  pendulum,  thus  they  can 
also  be  used  to  describe  the  motion  of  the  electrons  within  an  optical  wavelength, 
showing  the  evolution  of  electron  phase  and  phase  velocity  (analogous  to 
position  and  angular  velocity  for  the  classical  pendulum).  Since  phase  velocity  is 
related  to  the  energy,  phase  space  plots  give  an  informative  picture  of  the  energy 
exchange  between  electrons  and  the  optical  field. 

The  following  figures  show  the  phase  space  paths  for  20  sample 
electrons,  injected  with  the  same  initial  phase  velocity,  and  having  uniformly 
distributed  initial  phase.  The  FEL  phase  space  paths  can  be  found  from  the 
expression  [6] 

v2  =  v2  +  2  |a|  [sin  (^+^)- sin  (^o +  £>)],  (111.23) 

while  the  separatrix  is  given  by 

k2  =  2 1 a | [  1  +  sin ( £"+$>)]  i.e.,  v=  ±^2|a|  1  +  sin (^'+^)]/2,  (111.24) 
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which  defines  a  closed  contour  inside  of  which  are  the  closed  orbits,  and  outside 
of  which  are  the  open  orbits.  It  passes  through  the  unstable  fixed  points 
(£v)=(-;r/2,0)  and  (£  v)=  (3^/2, 0),  and  its  peak  to  peak  height  is  easily  seen 

to  be  4|a|1/: .  An  electron  moving  in  an  open  orbit  corresponds  to  a  pendulum  that 

swings  over  the  top,  whereas  an  electron  that  follows  a  closed  orbit  corresponds 
to  a  pendulum  that  swings  back  and  forth. 

In  Figure  7  it  is  plotted  the  evolution  of  20  sample  electrons  starting  at 
resonance  ( v0  =  0 ),  during  the  time  period  from  r  =  0  (appearing  as  yellow  dots) 
to  r=l  (appearing  as  red  dots).  All  paths  are  inside  the  separatrix,  thus  all 
electrons  follow  closed  paths.  The  symmetry  that  appears,  indicates  that  half  of 
the  electrons  gain  energy  from  the  optical  field  and  their  paths  shift  up,  while  the 
other  half  lose  energy.  This  is  represented  in  the  same  figure  from  the  graph  of 
the  net  gain  G(r),  which  is  almost  zero.  Gain  is  defined  as 

G(r>P^l~P°'  (IN.25) 

1  o 

where  7>(r)  is  the  optical  power  at  time  r  and  P0  is  the  initial  power.  The  optical 
phase  evolution  <p(r)  is  plotted  in  the  same  figure. 
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0.2 


Figure  7.  Phase  Space  evolution  of  Electrons  at  Resonance  ( vo  =  0), 
corresponding  optical  field  Gain  and  Phase  evolution. 

Contrary  to  the  nearly  zero  net  energy  transfer  when  electrons  are  injected 
at  resonance,  we  can  arrange  to  have  electrons  give  net  energy  to  the  optical 
field.  Injecting  electrons  slightly  off  resonance,  as  shown  in  Figure  8,  we  notice 
that  while  some  of  them  perform  open  orbits,  at  the  end  of  the  undulator  when 
t=  1,  most  electrons  have  moved  to  a  lower  value  of  phase  velocity,  resulting  in 
a  net  loss  of  energy.  In  this  case,  a  very  important  effect  takes  place  near  £=  n , 
namely  “bunching”  of  the  electrons  and  a  final  gain  of  more  than  10%. 

0.2 


o 

o.i 

o 


-k/2  5  3ic/2  0  t  1 

Figure  8.  Phase  Space  evolution  of  Electrons  injected  off  resonance 
(v0  =2.6),  corresponding  optical  field  Gain  and  Phase  evolution. 
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B.  DIFFRACTION  OF  LASER  BEAMS 


We  will  now  study  the  propagation  of  laser  light  through  a  sourceless 
medium,  in  order  to  be  able  to  understand  and  describe  how  light  travels  and 
moreover,  how  is  affected  by  diffraction.  We  start  with  Maxwell’s  equations, 
assuming  for  our  case  an  isotropic  and  homogeneous  medium.  When  there  are 
no  sources,  that  is  to  say  the  charge  density  p=  0  and  the  current  density  j  =  0 , 
we  write  the  four  Maxwell’s  equations  in  cgs  units  as 


V-E  =  0, 

(111.26) 

v-b  =  o, 

(111.27) 

w  p  1  SB 

VxE  = - , 

c  dt 

(111.28) 

Vx5  =  -  — , 
c  dt 

(111.29) 

where  c  is  the  speed  of  light.  Taking  the  curl  of  both  sides  of  equation  (111.28)  we 
have 

Vx|vx£ 


Applying  a  vector  identity  on  the  left-hand  side  and  substituting  VxB  from 
(111.29),  we  get 


/  — \  1  d 

'  1  dE ' 

(V-E  -V'£  = - 

1  ’  cdt 

c  dt  , 

(111.30) 


Then,  since  V -E  =  0  from  (111.26),  we  get 
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E  =  0, 


(111.31) 


f 

V 


V2 


c 2  dt2  j 


which  is  the  wave  equation  that  describes  light  propagation.  Electric  field  E  =  Es , 

where  a  is  the  polarization  vector,  is  a  function  of  position  and  time,  thus 
equation  (111.31)  implies  four  second-order  derivatives  in  (x,t).  Due  to  the 

coherence  of  the  laser  light,  electric  field  E  is  slowly  varying  in  the  direction  of 
propagation  over  an  optical  wavelength.  Likewise,  E  is  slowly  varying  in  time 
over  the  optical  period.  Then,  we  can  write  the  magnitude  of  the  complex  electric 
field  in  the  form 

E  (x,r)  = 

=  |a(jc,f)|e,*sVfe““')I  (111.32) 


where  kz-cot  is  the  carrier  wave,  a  =  \a\eip  is  the  complex  optical  field  and 
x  =  (x,y,z).  Next,  we  want  to  substitute  the  electric  field  (111.32)  into  the  wave 
equation  (111.31).  After  taking  the  required  derivatives,  we  find 


V  ^GC  + 


( 82a  da  ,  2  'I 

— -  +  2ik - k  a  - 

V  dz2  5z 


J-7l 


1  (  d~a 


'  da  2  ^ 

— —-hat - a>  a 

dt2  dt 


i(kz—a)t )  _ 


0, 


1.33) 


with  V2  =d2/dx2  +d2/dy 2  the  transverse  Laplacian.  We  can  simplify  this  equation 
by  dropping  out  the  second  derivatives,  because  laser  field  amplitude  |a(x,r)| 
and  phase  cp{x,t)  are  slowly  varying.  This  means  that 


d  a 
dz 


=  k\a\  , 


d  a 
dt 


co  a 


S(P  j 

—  =  cocp  and 

dt 


(111.34) 
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From  d\a\/dz  =  k\a\  we  find  that  d2\a\/dz2  =  kd\a\/dz  and  similarly  from 

8\a\/dt=  co\a\  we  find  that  d2\a\/dt1=  cod\a\/dt,  so  that  we  are  justified  in 

dropping  the  second  derivatives  in  favor  of  the  first  derivatives  in  (111.33).  In  the 
same  equation,  the  last  term  in  each  parenthesis  cancels  out,  since 

k=  co / c=>£2  =  co2  / c2  =^>  k2a  =  co2a / c2 .  Then,  after  multiplying  by 
equation  (111.33)  becomes 


( d 

i  dYI 

Vl+2  ik 

— +  — 

{dz 

c  dt  J 

a(x,t )  =  0 . 


(111.35) 


In  order  to  simplify  the  last  expression,  we  introduce  a  new  variable 
u  =  z-ct  and  again  use  the  dimensionless  time  z  =  ct/L  (where  L  is  the 

propagation  range).  We  can  now  rewrite  the  operator  d/dz  +  c~l  d/dt  that  appears 
in  (111.35)  taking  the  necessary  partial  derivatives 

d  _du  d  dz  d 

dz  dz  du  dz  dz  >  (111.36) 

1  d  _  1  du  d  1  dz  d 
cdt  c  dt  du  c  dt  dz 


and  since  du/dz  =  1,  dz/dz  =  0,  du/dt  =  -c  and  dzldt  =  c/L,  the  operator 
d/dz  +  c^d/dt  is  just  Uld/dz  and  consequently  the  wave  equation  (111.35)  is 
now  written 


f 


vl  + 


v 


2 ik_d_' 
L  dz j 


a(x,t)  =  0 , 


(111.37) 


which  is  called  the  “parabolic”  or  “paraxial”  wave  equation,  because  we  have 
made  use  of  the  paraxial  approximations  in  (111.34). 
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At  this  point,  it  is  convenient  to  use  dimensionless  transverse 
coordinates  (x,y)  along  with  the  normalized  time  z .  If  we  rearrange  the  terms  in 
equation  (111.37),  we  can  rewrite  it  as 


r  -iC 

l  ~2ky 


vi  + 


A 

8z 


a(x,r)  =  0 , 


(111.38) 


where  we  notice  from  the  first  term,  which  is  the  diffraction  term  of  the  wave 
equation,  that  the  transverse  extent  of  the  laser  beam,  together  with  the 
wavenumber  k  =  27z/A,  and  the  propagation  range  L,  affect  the  diffraction  of  the 
beam.  The  equation  takes  a  simpler  form  if  we  normalize  the  transverse 
coordinates  using  the  relationships 


x  =  x 


f  k  > 

Yi 

f  k  \ 

v2  L) 

j 

y  =  y\ 

(111.39) 


In  (111.39),  the  normalization  factor  (A:/2l)^  is  related  to  the  characteristic  mode 
waste  radius  WQ  by 


W„  = 


(111.40) 


with  Za  =  L  in  our  case,  where  Zc  is  the  Rayleigh  length,  defined  to  be  the 
propagation  distance  in  which  the  beam  doubles  its  area  and  it  is  determined  by 
the  curvature  of  the  mirrors.  Finally,  using  only  dimensionless  coordinates  and 
utilizing  the  dimensionless  transverse  Laplacian,  we  arrive  at  the  dimensionless 
parabolic  wave  equation,  written  without  the  tildes  for  convenience 


f 

v 


JA 

dr  j 


a(x,r)  =  0 . 


(111.41) 
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From  this  last  wave  equation,  we  notice  that  for  a  beam  with  transverse 
dimensions  of  the  order  of  unity  or  less,  we  have  to  take  into  account  the 
diffraction,  since  it  affects  the  beam  a  great  deal,  over  the  range  r=  0-»  1 ,  while 
for  a  much  wider  beam  we  can  ignore  it.  Diffraction  of  the  wavefront  causes  both 
the  amplitude  and  the  phase  of  the  optical  field  to  evolve. 


C.  THE  FEL  WAVE  EQUATION 


We  have  found,  so  far,  a  mathematical  description  for  the  laser  light  that 
propagates  in  absence  of  sources  and  also  for  the  electron  motion  in  the 
presence  of  light.  What  is  left,  in  order  to  generalize,  is  to  describe  the  optical 
wave  evolution,  taking  into  account  the  interaction  with  the  current  source  due  to 
the  electron  beam.  We  start  with  the  full  wave  equation  [6] 


f 

v 


V2 


c2  dt2  j 


A[x,t ) 


4n-  , 

- —  J±(X,t)  , 


(111.42) 


where  A(x,t)  is  the  optical  vector  potential  and  J±(x,t)  is  the  transverse  current 

density  due  to  oscillations  of  the  electrons  passing  through  the  undulator.  Then 
we  are  able  to  find  the  optical  electric  and  magnetic  fields  from  the  vector 
potential  using  the  relationships 

—  1  d~A  —  — 

£  =  B=VxA.  (111.43) 

c  dt 

Once  again,  we  will  make  use  of  the  fact  that  laser  light  is  coherent.  Thus,  we 
take  A(x,t )  to  have  a  slowly  varying  envelope  in  the  direction  of  propagation.  In 

addition,  we  take  the  electric  field  to  vary  slowly  over  a  laser  wavelength  and  to 
vary  slowly  in  time,  relative  to  the  optical  frequency.  Assuming  a  complex  laser 
electric  field  E  =  \E\el,p ,  the  vector  potential  can  be  written  [6] 
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(111.44) 


A(x,t ) 


E  (x,t) 


s  e 


i{kz—cot ) 


where  kz-  cot  is  the  carrier  wave  and  a  is  the  laser  polarization  vector,  which  is 
s  =  (—i,l,0)  for  the  case  of  a  helical  undulator  producing  circular  polarized  laser 

light.  If  we  assume  again  slowly  varying  wave  amplitude  and  phase  and  use 
(111.44),  the  wave  equation  (111.42)  becomes 


g  ^i(kz-cot) 


k 


( a 

i  a\] 

V3 +2  ik 

—+ — 

ydZ 

c  dt  J 

c 


_ _ _  o  _  ,  1 

( d 

i  a  VI 

V3  +27*1 

— + —  1 

l ydz 

c  dt  J  \ 

—i 


e 


\kz-cot) 


(111.45) 


We  will  now  use  the  same  technique  as  we  used  in  (111.36),  introducing  the  new 
variable  u  =  z-ct  and  r  =  ct/L  ,  where  L  is  the  length  of  the  undulator,  so  as  to 

rewrite  the  operator  d/dz  +  c  1  d/dt .  Then  (111.45)  becomes 


vi  + 


2  ik  d 
L  dr 


e=-™Lj1 


s  e 


i{kz-cot ) 


(111.46) 


This  last  equation  is  the  parabolic  wave  equation  with  a  source  current 
J±,  due  to  the  bunched  electrons.  The  source  current,  which  is  due  to  all  of  the 
electrons,  can  be  expressed  as 

(111.47) 


where  ^3)  is  the  three-dimensional  Dirac  delta  function,  defined  as 
S(3) (r)  =  S(x)S(y)S(z)  and  t )  is  the  position  of  the  ith  electron  at  time  t.  The 

transverse  velocity  /?,  found  in  (111.10)  can  be  expressed  as 
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(111.48) 


( cos kQz, sin kaz, 0)  =  5Kei -  —  ise  ,k°z 

r  {  r 


Plugging  (111.48)  into  the  source  current  relationship  in  (111.47)  and  then  the 
resulting  source  current  into  the  parabolic  wave  equation  (111.46),  we  have 


2 ik  d  ' 
+- 


L  8z 


E  =  -4 nieKk  p(x,t )  (< ei( /y )  (J  f) . 


(111.49) 


In  this  equation,  p(x,t)  is  the  electron  beam  particle  density  in  a  small  volume 
element  dV 


P  (*>  0  =  J  X  ^(3)  “  A  (0)  dV , 

V  i 


(111.50) 


and  with  the  symbol  ( . )  we  denote  the  average  over  sample  electrons. 

Employing  again  the  dimensionless  time  r  =  cr/7,  where  L  is  the  length  of  the 

undulator,  and  defining  the  dimensionless  complex  laser  field  a  =  |a|<?'?’,  with  the 
laser  field  amplitude  |a|  =  4 xNeKLE / y20mc2 ,  as  used  in  (111.22),  the  wave  equation 
becomes 


f 

v 


2k  dz 


i(x,r)  =  - 


(111.51) 


where  we  have  introduced  another  important  FEL  parameter,  the  dimensionless 
current  density  j  =  87r2Ne2K2L2p/y3omc2 ,  and  we  have  made  the  reasonable 
assumption  that  during  the  whole  interaction.  We  want  now  to  make  the 

wave  equation  completely  dimensionless,  thus  we  bring  in  again  the 
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dimensionless  transverse  coordinates  x  =  x(£/2Z)^,  y  =  y{kj2L')//'1  and  the 

dimensionless  Laplacian  operator.  Dropping  the  tildes,  the  FEL  wave  equation  is 
finally  written 


f 

V 


A 


— V:  +  — 
4  dr 


1.52) 


(x,r)  =  - 


In  the  case  when  the  diffraction  can  be  neglected,  the  last  equation  is 
simplified  to 


da  _°  _ 
dr 


(111.53) 


From  this  expression,  we  realize  that  the  electron  beam  should  be  bunched, 
since  the  factor  (e~K)  is  a  measure  of  the  electron  bunching.  For  instance,  if  the 
electrons  are  bunched  at  around  £=n,  the  laser  field  amplitude  goes  up, 
whereas  if  the  bunching  occurs  at  £=  0 ,  the  field  amplitude  goes  down.  Equation 
(1111.53)  also  suggests  that  the  laser  field  is  related  to  the  dimensionless  current 
density  j ,  so  that  when  j  <  n  we  are  in  the  low  gain  regime,  while  for  j  »  n  the 
coupling  between  the  laser  light  and  the  electron  beam  is  large  and  we  have  high 
gain. 


D.  GAIN 

In  order  to  better  understand  and  describe  the  energy  exchange  between 
the  electron  and  the  optical  beam,  it  is  handy  to  make  use  of  the  gain  we  defined 
in  (111.25),  which  is  the  fractional  change  of  the  optical  field  power  during  a  single 
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pass.  In  terms  of  the  optical  field  amplitude,  gain  can  be  written  as 


G(r)  = 


a<xf-v  K*or 


-l 


or 


G(t)«  2 


i(r)|  ^ 

—  -II, 


(111.54) 


for  low  gain  and  small  change  in  |a(r)|  from  the  initial  field  amplitude  a0. 

It  is  particularly  useful  to  study  the  gain  in  the  weak  optical  field  state, 
because  it  gives  us  information  about  the  performance  of  the  FEL.  For  that 
reason,  we  will  derive  an  expression  for  the  gain,  assuming  weak  optical  fields, 
a a=  7i  and  low  gain,  so  as  to  take  |a|«a0and  0  for  the  whole  interaction. 
Applying  the  perturbation  theory,  we  start  by  expanding  the  electron  phase  C, 
and  the  electron  phase  velocity  v  in  powers  of  the  initial  small  field  amplitude  a0 

^=^0)  +  ^1)  +  ^2)  +  ...;  1=  v=  i/°>  +  v(1)  +  v(2)  +  --->  (111.55) 

where  is  the  zeroth  order  in  aa,  ^  is  the  first  order,  and  so  on.  We  can  now 
plug  (111.55)  into  the  pendulum  equation  (111.22)  to  get 

<£'(°°)+<^(°1)+  <^(24 —  =  v(0)+  v^+  v,2*h —  =  ao  cos^°*  +  ^(I)  h — j .  (111.56) 


From  the  zeroth  order  in  a0  we  have  ^(0)=k(0)=0  which  gives  the  solutions 

<^0)  =  +  v0r  =  0  and  v(0)  =  vo .  From  the  first  order  in  ao ,  using  the  zeroth  order 

solutions  and  integrating  overtime  r,  we  find  the  solutions 


=  -^r[c°s(C  +  v0t)~  cos(c)+  v0rsm(C0)], 

^ o 

•/0  =  f  =  -[sin(C  +  v,t)- sin(C)] .  (111.57) 

^ o 
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We  can  now  make  use  of  these  results  in  the  similarly  expanded  wave  equation 
(111.53) 

a  =  -  j(e~i()=  -  j(e-i^)+^2)+  -)) .  (111.58) 

Taking  into  account  only  the  zeroth  and  the  first  order,  we  can  approximate  the 
last  equation 

a  *  -  j{e-^e-^)K  -  j{e^\  1  -  i^)) ,  (111.59) 

but  since  (e~i^">)  =  {e~ii(o+v°T))  =  J  e~i(io+v°T)d^0  jin  =  0 ,  (111.59)  is  written 

a  «  jie-^uZ®)  =  ij{eii<:°+v°T)f')) .  (111.60) 


We  can  plug  in  ^  from  (111.57)  and  integrate,  recalling  that  the  complex  optical 
field  is  a  =  |a|e'v,  to  find 


la (r)|  =  all  +  /  2  ~  2 c°s(nO-  v.rrinM])  + . . . , 


¥s)  =  j 


H1+4 - 

(  2sin(vor)-  v0r( 1  +  cos(vor)Jj 


Now,  using  equation  (111.54),  we  are  able  to  arrive  at  the  low  gain  equation  along 
the  undulator  [6] 


G(r)  =  j 


f  2  - 2cos(t/0r)-  v0rsin(v0r)') 

l  v03  y 


(111.62) 
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This  equation  indicates  that  the  gain,  in  the  weak  field  approximation,  is 
proportional  to  the  dimensionless  current  j  and  the  initial  electron  phase  velocity 
vo .  Plotting  the  final  gain  spectrum  at  r  =  1  in  a  single  pass  through  the 
undulator,  as  shown  in  Figure  9,  we  can  verify  the  results  we  found  in  the  phase 
space  plots  of  Figures  7  and  8.  Particularly,  when  the  electrons  are  ejected 
exactly  at  resonance  vo  =  0,  then  there  is  zero  net  gain,  meaning  that  the 
electrons  lose  and  gain  the  same  amount  of  energy  during  their  interaction  with 
the  optical  field  inside  the  undulator.  We  also  note,  that  the  gain  curve  is  anti¬ 
symmetric  in  vo .  The  peak  gain  of  about  13%  is  obtained  when  the  FEL  operates 
off  resonance  at  v0  ~  2.6 ,  while  the  peak  absorption  occurs  at  vo  «  -2.6  .  These 
results  are  consistent  with  the  phase  space  analysis,  suggesting  that  we  must 
operate  the  FEL  slightly  off  resonance  in  order  to  have  the  best  exchange  of 
energy  from  the  electron  beam  to  the  optical  mode  and  hence  the  maximum 
gain. 


Figure  9.  Single  pass  gain  versus  initial  electron  phase  velocity  va  at  r  =  1. 
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IV.  OPTICAL  THEORY 


Up  to  this  point  we  have  examined  the  interaction  of  the  electron  beam  and 
the  optical  field  inside  the  undulator,  arriving  at  some  important  expressions.  We 
have  also  introduced  the  gain  G(r) ,  a  useful  tool  to  measure  the  performance  of 
this  interaction.  In  this  section,  we  will  focus  our  study  only  on  the  optical  field, 
propagating  in  free  space.  Equation  (111.52)  describes  the  optical  field  evolution, 
in  the  more  general  case  of  the  presence  of  source.  However,  when  light  is 
traveling  in  a  sourceless  medium,  specifically  inside  the  undulator  while  there  is 
no  electron  beam  present  and  outside  the  undulator,  this  equation  reduces  to 
equation 


f 

v 


8  ' 
dr  j 


a 


(x,r)  =  0 . 


(IV.1) 


In  the  following  sections,  we  will  find  solutions  to  this  wave  equation,  starting  with 
the  one  that  describes  typical  optical  waves  produced  by  a  Free  Electron  Laser. 
We  call  this  solution  the  “fundamental”,  and  it  is  the  lowest-order  mode.  Higher- 
order  modes,  which  we  are  going  to  analyze  later,  may  also  appear. 


A.  FUNDAMENTAL  SOLUTION  TO  THE  WAVE  EQUATION 

An  ideal,  realizable,  laser  beam  cross  sectional  shape  is  Gaussian.  We 
therefore  seek  a  Gaussian  solution  to  the  wave  equation.  A  trial  solution  of  this 
form  can  be  written 


a(r,r)  =  aop(r)e  9<r) ,  (IV. 2) 

where  a (r,r)  is  the  dimensionless  optical  field,  and  aG  is  the  optical  field 
amplitude  at  the  point  (x,y,r)=  (0,0,0),  r2  =  x2  +  v2 .  The  functions  p{ r)  and  q( r), 
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as  yet  unknown,  will  be  determined  by  the  wave  equation  (IV. 1)  and  the  initial 
conditions  imposed  by  the  desired  Gaussian  shape  of  the  beam, 


a(F,0)  =  aoe  z° , 


(IV  .3) 


where  zo  =  ZJL  is  the  dimensionless  Rayleigh  length,  and  L  is  the  propagation 
distance.  These  initial  conditions  give 

/?(0)  =  1  and  q(0)  =  zo.  (IV. 4) 


In  order  to  find  the  expression  we  are  looking  for,  from  the  trial  solution  (IV. 2),  we 
have  to  find  the  functions  p(r)  and  q{ r).  Before  putting  (IV. 2)  into  the  wave 
equation,  we  switch  to  cylindrical  coordinates,  due  to  the  symmetry  of  our 
problem.  Then,  the  wave  equation  (IV. 1)  becomes 


i  1  d 

r  d N 

4  r  dr 

V 

a 

+ — 
dr 


a(r,r)  =  0  . 


(IV  .5) 


Now  we  insert  the  trial  solution  (IV. 2)  into  (IV. 5).  The  first  derivative  of  a (r,z) 
with  respect  to  r  is 


«i  =  _2a  ^,*1, 

dr  q(r) 


(IV  .6) 


from  which 


If  3]  =  4a 


dr\  dr)  0  q(r)  \  q(T) 


'  'J  -.1. 


(IV.7) 


Finally,  we  calculate  the  time  derivative  of  a (r,z) 
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da 


=  a0p(r)e 


<?(r) 


f  1  dp  r  dq\ 


■  + 


dz  U1  '  '  l  p(z)  dz  q\z)  dz 


(IV  .8) 


We  can  now  go  back  to  the  wave  equation  (IV.5)  and  employ  these  derivatives 


i  1  d  f  da 

- r —  +  —  =  0 

4  r  dr  v  dr)  dz 


.  pif)  ~dd)(  r2  7)(  1  dP  ,  }'2  dci\ 
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°  q(z)  \q(r) 
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-1 1  +aop(z)e 


■  +  ■ 


p(z )  dz  q  (z)  dz. 


=  0 


(IV  .9) 


Dividing  through  by  a 0p(z)e  /q(T)  we  get 


r  1  1  dp  r  dq 

-l -  +  1 - + - —  +  — : - -  =  0 


q(z)  q(z)  p(z)  dz  q\z)  dz 

Multiplying  by  -iq2(z)  and  rearranging,  we  can  write  the  last  expression  as 


(IV.10) 


2i  ,  .dq  |  ,  2,  s 
r  \-l-i  —  \+q(T) 

dz 


I  i  dp ^ 
\q(z)  p(z)dz) 


=  0 


(IV.11) 


In  order  for  this  equation  to  hold  for  all  r  at  any  time  z  ,  the  quantities  in  the  two 
parentheses  must  be  identically  zero,  and  therefore 


.  dq 
dz 


1  i  dp 


q(z)  p(z)  dz 


and 

(IV. 12) 

0  . 

(IV. 13) 

From  equation  (IV. 12),  after  integrating,  we  find  q(z)  =  iz  +  q0  ,  and  from  the  initial 
conditions  (IV. 4)  we  see  that  q(0)  =  zo,  so  that 
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q(T)  =  ir  +  zo  . 


(IV. 14) 


Using  the  last  result  in  equation  (IV.  13)  we  see  that 

dp  _  ip(r)  ^  p(r) 
dz  iz  +  za  -r  +  iz0 

\—dp  =  \^—=. 

*>(0  ]~T'+iZ0 

In p( t)  =  -  ln(- z '+  izo )  |q  =  - ln(- r  +  izo )  +  ln(/zo )  ,  (IV.  15) 


and  finally 


p(r)  = 


1 +  1T/ 


(IV. 16) 


We  can  now  substitute  p( r)  and  q(r)  into  (IV.2)  and  write  down  the  solution  to 
the  parabolic  wave  equation  as 


i(/%r)  =  ; 


\+iV 

V  /ZoJ 


(IV. 17) 


This  solution  can  be  written  in  a  more  informative  way,  expressed  in  terms  of  the 
dimensionless  beam  width  w(z),  which  is  defined  as 


( .  z1  \  N  |  z~ 


w(*)  =  Azo\\  +  ~ rj  or  w^  =  yjzo  +  Y- 


(IV. 18) 


We  note  that  at  r  =  0 ,  the  beam  width  or  waist  in  this  case  (the  minimum  value  of 
the  beam  width),  is  just  w(0)  =  wo  =  Jz^  .  Of  course,  the  beam  waist  need  not  be 
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at  the  beginning  of  propagation  distance,  thus  a  more  general  expression  for  the 
beam  width  at  any  point  is 


w(r)=r+^^,  (iv. 19) 

where  rw  is  the  dimensionless  position  of  the  beam  waist  (normalized,  as  usual, 
to  the  propagation  distance  L  ).  A  cross-section  of  a  propagating  Gaussian  beam 
is  shown  in  Figure  10.  Physically,  the  beam  width  is  the  distance  where  the 
optical  field  drops  to  1/e  of  its  maximum,  as  shown  in  Figure  1 1 . 


Figure  10.  Cross-section  of  a  propagating  Gaussian  Beam. 

We  are  now  ready  to  go  back  and  modify  the  solution  (IV.  17).  To  make  the 
calculations  easier  we  let  u  =  \  +  ixjzo .  Then,  equation  (IV.  17)  becomes 


a (r,z) 


u 


(IV. 20) 
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Figure  1 1 .  Gaussian  beam  profile. 


and  since  the  complex  conjugate  of  u  is  u*  =  1  -  ir/z0  ,  we  can  write  (IV. 20)  as 


i(/“,r)  =  ao  —e 


a(r,r)  =  ao-e 
u 


r2  (\—iz/z0) 


r  izr' 


'ho 


n  ( —  \  ^  1  znuu  znuu 

a(r,r)  =  ao-e  e 
u 


(IV. 21 ) 


The  complex  number  u  can  be  written  in  polar  form  as 


u  = 


/  .  r t') 

I  T  )  i  arctanl  - 

1+id e  ■ 


(IV. 22) 


and  since  uu*  -  l  +  (r/z0)2  ,  equation  (IV.21)  becomes 
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a(/_:,r)  =  ao 
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e  x 


1  +  (r/z<J 


I-  l  TV 

fc)^.('-!/v)  -v(i+r/v) 


(IV. 23) 


We  now  insert  the  beam  width  w(t),  defined  in  (IV.  18),  into  the  last  equation, 
which  becomes 


i(r,r)  =  . 


w„ 


w  O) 


-/arctan  ^ 


°  Mr) 


(IV. 24) 


where  we  also  used  wa  =  y[z~0  .  Finally,  we  can  write  the  solution  to  the  parabolic 

wave  equation  (IV.  1)  more  compactly,  realizing  that  the  last  two  exponents  in 
(IV. 24)  are  pure  imaginary.  So,  we  let 

_  2  f 


t)  = 


rr 


2  ,  2 

Z+T 


-  arctan 


\ZoJ 


and  write  the  fundamental  solution  as 


(IV. 25) 


i(/%r)  =  ; 


w„ 


w(t) 


V(f) 


(IV. 26) 


B.  PROPERTIES  AND  PROPAGATION  OF  GAUSSIAN  BEAMS 

The  Gaussian  beam  is  the  lowest-order  mode.  A  beam  of  this  profile  can 
be  fully  described  by  the  beam  waist  w0  and  the  position  of  the  beam  waist  rw . 
We  will  study  the  behavior  of  the  optical  beams  described  by  the  solution  (IV.26), 
in  order  to  better  understand  the  evolution  of  their  properties  as  they  propagate. 
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First,  we  note  that,  at  the  beginning  of  the  propagation,  when  z  =  0,  the 
beam  has  the  Gaussian  shape  that  we  chose  with  the  initial  conditions,  namely 

a(r,0)  =  aoe  z° .  (IV.27) 

As  time  increases,  the  beam  width  defined  in  (IV.  18)  gets  bigger  and  so 
the  transverse  area  of  the  beam  expands.  As  a  consequence,  for  a  fixed  distance 
from  the  center  (fixed  r ),  the  field  strength  goes  down  as  the  optical  field 
propagates.  Looking  at  the  center  of  the  optical  mode,  where  r  =  0 ,  the  field 
decreases  continuously,  as  it  is  expected,  since  a(0,r)  goes  like  ~  1/r . 

Gaussian  modes  have  the  very  important  feature  that  they  do  not  change 
their  profile  as  they  propagate,  except  that  the  beam  radius  varies.  In  other 
words,  a  Gaussian  beam  stays  Gaussian.  This  is  the  reason  why  we  where 
looking  for  a  solution  to  the  wave  equation  of  this  form. 

Figure  12  represents  a  computer  simulation  of  the  propagation  of  a  laser 
light  in  the  fundamental  mode.  In  this  figure  is  plotted  a  cross-section  of  the 
Gaussian  mode  as  it  propagates  from  r  =  0  to  z  =  1 ,  together  with  the  transverse 
cross-section  at  the  beginning  and  at  the  end  of  the  propagation  distance. 


Figure  12.  Cross-section  of  a  propagating  Gaussian  Beam  with  zw  =  0  and 

zo  =  0.3 . 

Additionally,  in  Figure  13  is  shown  the  surface  plots  of  a  Gaussian  mode,  at 
the  beginning  of  the  propagation  when  z=  0,  at  r  =  0.5  and  at  the  end  of  the 
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propagation  distance  when  t=  1  (note  the  difference  in  the  vertical  scale 
between  the  three  graphs).  From  all  these  plots,  we  are  able  to  visualize  what  we 
analyzed  previously,  specifically  the  fact  that  the  light  beam  expands  as  it 
propagates  and  also  that  it  suffers  attenuation. 


(a) 


(b) 


Figure  13.  Surface  plots  of  a  Gaussian  Beam  evolution  at  times  (a)  r  =  0, 
(b)  t  =  0.5  and  (c)  t  =  1 .  Scales  in  z-axis  are  different. 
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V.  HIGHER  ORDER  MODES 


In  the  previous  chapter,  we  considered  only  the  fundamental  Gaussian 
solution  to  the  wave  equation,  which  has  the  simplest  intensity  profile  and  is  a 
good  approximation  to  describe  a  laser  beam.  However,  if  we  want  to  generalize 
for  non-ideal  beams  and  include  the  possible  variations  in  the  transverse  beam 
profile,  we  have  to  find  other  solutions  that  correspond  to  higher  order  modes. 
These  solutions  to  the  parabolic  wave  equation  have  similar  properties  to  the 
Gaussian.  Furthermore,  they  form  a  complete  and  orthogonal  set  of  functions, 
called  the  modes  of  propagation,  so  that  every  arbitrary  optical  beam  can  be 
expanded  in  terms  of  them  [7],  We  will  find  that  the  fundamental  Gaussian  beam 
is  just  the  lowest-order  mode  of  the  generalized  solution. 

As  one  can  guess,  the  need  for  studying  the  optical  beam  more 
extensively  arises  from  both  experimental  results  and  simulations.  Generally 
speaking,  whenever  an  FEL  operates  under  non-ideal  conditions,  the  laser  beam 
may  contain  more  than  one  higher-order  modes.  Such  conditions  could  be 
vibrations  of  the  mirrors  in  a  portable  system,  electron  beam  shifts  and  tilts, 
mirror  imperfections  or  deformations,  and  so  on.  For  instance,  vibrations  can 
cause  the  electron  beam  to  be  injected  slightly  off-axis  and  to  deform  the  mirrors, 
while  heat  can  change  the  radius  of  curvature  and  shape  of  the  mirrors. 


A.  HERMITE-GAUSSIAN  MODES 

The  Hermite-Gaussian  solutions  are  derived  from  the  parabolic  wave 
equation  in  rectangular  coordinates,  starting  with  the  more  general  trial  solution 
[8] 


a(r,z-)  =  a0  g 


(  „  A 


\w(t)  j 


f  y  A 


P(  !■)+' 


z„  9(0 


(V.1) 
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where  the  functions  p{r),  q{r),  g(x/w)  and  h(y/w)  are  found  following  a 

procedure  similar  to  the  one  used  for  the  fundamental  mode.  The  same  initial 
conditions  as  before  require  that 

/>(0)  =  0,  <?(0)  =  1,  (V.2) 

g(x/w0)=l,  h(v/w0)=  1. 

Using  (V.2)  and  substituting  the  trial  solution  into  the  wave  equation,  we  arrive  at 
equations  that  have  the  form  of  the  Hermite  equation  and  can  be  solved  in  terms 
of  the  Hermite  polynomials.  After  finding  all  the  four  unknown  functions  p,  q,  g 
and  h  ,  we  write  the  m,n  mode  of  the  Hermite-Gaussian  beam  as 


(A'' 

H 

(Ay) 

{  w  J 

n 

{  w  J 

w 

(V.3) 

where  the  phase  is 

l'^T 

0  ,  -  ( m  +  n  +  I ) arctan 

w~w~ 


e  A v2  et™Ar,T) 


f  T  ' 


\ZoJ 


(V.  4) 


and  m,n  are  positive  integers  identifying  the  mode.  In  (V.3),  Hmand  H„are  the 
Hermite  polynomials  with  m,  n  the  non-negative  integer  indices.  These  indices 
determine  the  shape  of  the  beam  profile  in  the  x  and  y  direction  respectively.  The 
first  ten  of  the  Hermite  polynomials  are  listed  in  Table  1.  It  is  obvious  that  when 
m  =  0  and  n  =  0  ,  the  Hermite-Gaussian  beam  is  just  the  Gaussian  beam,  which 
is  why  it  is  called  the  fundamental  mode.  A  few  examples  of  higher  order 
Hermite-Gaussian  modes  are  shown  in  Figure  14,  where  we  note  that  m  nodes 
appear  running  vertically,  while  n  nodes  run  horizontally. 
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tf„60=i 

Hx  (m)=  2u 
H2(u)=  4u2-2 
//3  (u)=  8u3  -  12m 
//4(m)=16u4-48m2  +  12 
(m)=  32u5  -  160m3  +  120m 
H6  (u)=  64u6  -  480m4  +  720 u2  -  120 
H1  (u)=  128u7  -  1344m5  +  3360m3  -  1680m 
Ha  (m)=  256u8  -  3584m6  +  13440m4  -  13440m2  +  1680 
(m)=  512u9  -  9216m7  +  48384m5  -  80640m3  +  30240m 


Table  1 .  Hermite  polynomials 


m  =  2,  77  =  0  777  =  1,  77  =  2  777  =  2,  77  =  3 

Figure  14.  Hermite-Gaussian  modes  in  the  transverse  plane  x,y.  Red 
indicates  highest  intensity,  while  blue  indicates  least  intensity. 
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B.  LAGUERRE-GAUSSIAN  MODES 


The  optical  cavity  of  an  FEL,  utilizing  circular  mirrors  and  a  round,  aligned 
electron  beam,  has  cylindrical  symmetry.  This  fact  suggests  that  we  should  use 
cylindrical  coordinates  to  solve  the  dimensionless  parabolic  wave  equation  (IV. 1), 
which  is  now  written 


d~a  1  da  1  d2a  A  da 

— 7  + - +  — 7  +  4/ 

dr2  r  dr  r~  d9~  dr 


(V.5) 


since  the  dimensionless  transverse  Laplacian  in  cylindrical  coordinates  is 


d2  1  d  Id2 

- 7  H - 1 - ~ - 7  . 

dr~  r  dr  r~  d& " 


(V.6) 


Motivated  from  the  success  in  finding  the  Hermite-Gaussian  modes  trying  the 
solution  (V.1),  we  will  now  put  a  similar  exponential  term  in  the  trial  solution  for 

,7  Q 

the  wave  equation,  adding  e  to  match  the  single  valued  periodic  boundary 
conditions,  where  /  is  an  integer.  We  also  guess  a  power  series  in  r,  and  since 
we  anticipate  diffraction  that  suggests  the  use  of  r/w  dependence,  we  try  the 
solution 


a(F,<9,r)  =  ac 


V2  r 

Mr) 
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2  r 


y«>. 


A  -  P{T)+^—-iW 


(V.7) 


where  t{i r'/w1)  ,  P(r)  and  q(r)  are  the  unknown  functions  that  we  will  try  to 
find  next.  Although  we  use  cylindrical  coordinates,  we  still  require  that  at  r  =  0 


a(r,i9,0)  =  aoe  z°e'w , 


(V.8) 
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which  gives  £{2/ /w2)  =  \,  P{ 0)  =  0  and  g(0)  =  l.  In  order  to  put  the  trial  solution 
(V.7)  into  the  wave  equation  (V.5),  we  first  need  to  calculate  the  derivatives  of  the 
function  £(2 P/w2)  using  the  chain  rule 


s((2r '/».')  _  at(  >■/»■)  a(2r!/».2)  _  4,. 

*  a(2r/wJ)  ar  w2  '  1  ' 

where  we  used  the  prime  to  represent  the  first  derivative  of  £  with  respect  to  its 
argument  2  r2 /w2  .  Similarly,  using  the  double  prime  for  the  second  derivative, 
we  find 


d2£(2  r2/w2) 

dr2 


d_(4 r 
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8£') 

+  r - 

dr  J 


The  time  derivative  of  £(2r2/w2)  is 


(V.10) 
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with  (..)  indicating  the  derivative  with  respect  to  the  dimensionless  time  z .  Now 
we  are  ready  to  calculate  the  derivatives  of  the  trial  solution  (V.7) 
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(V.14) 
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(V  .16) 


Next,  we  want  to  substitute  the  trial  solution  (V.7)  into  the  wave  equation  (V.5). 

/  i—  i  \i  \p+—- -<■«] 

Bringing  all  the  derivatives  together  and  dividing  through  by  a0U2/w\ey  z°q  J 
we  have 
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r  and  i  dependence 


only  dependence  no  r  dependence 


We  want  this  equation  to  have  solutions  for  each  transverse  plane,  or  in  other 
words,  for  each  fixed  time  z .  Thus,  we  need  to  set  the  r  and  £  dependent  terms 
equal  to  a  constant,  while  the  r 2  dependent  terms  should  be  equal  to  zero 
because  equation  (V.17)  must  be  valid  for  any  transverse  position.  Eventually, 
the  terms  with  no  r  dependence  must  be  equal  and  opposite  to  the  constant 
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value  of  the  terms  dependent  on  r  and  £ ,  in  order  for  the  equality  to  be  true. 
Hence,  we  have  to  solve  the  following  three  equations  simultaneously 


4rf  2l_t  4r  f+2f  4  iwr2  £' 
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(V.18) 

(V.19) 

(V.20) 


where  we  used  -  £  as  the  constant  for  the  r  and  £  dependent  terms.  We  will 
start  from  (V.19),  since  it  is  straightforward 

2  •  2  °  2  ( 

r  ir  q  n  r  I 

— — 7  +  - 2  =  0  4>  - ^-1 

4  q  q"  za  q-  K 

where  for  the  last  step  we  use  the  same  argument  as  we  did  for  the  derivation  of 
the  fundamental  mode,  namely  that  the  equation  must  hold  for  all  r.  After 
integrating,  and  since  #(0)  =  1  from  the  initial  conditions,  we  find 


1  °^i  1  0 

- f  i  q  I  =  0  => - f  i  q  =  0  ,  (V.21 ) 

z„  )  z„ 


q(T)  =  -  +  1.  (V.  22) 
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From  the  definition  of  the  beam  width  w(t)  =  s]zo  +  t2!zo  in  (IV.  18),  we  have 
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w  =  t / wzq  .  Rearranging  terms  in  equation  (V.18)  we  write 
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where  the  terms  in  the  parenthesis,  after  substituting  q(r)  from  (V.22)  and  using 
z  =  w2 ,  become 


2 r2  zo  q  2 r2  zo  2 r2  2 r2 

Let  u  =  2r2 jw2  to  rewrite  (V.23)  more  compactly  as 


(V.24) 


n\ "  +  (/  +  1-m)I  '  +  ~^\  =0,  (V.25) 

which  has  the  form  xy"  +  (l  +  1-  x)y'  +  py  =  0  ,  of  the  differential  equation  for  the 
associated  Laguerre  function  [9],  where  /  and  p  are  real  numbers.  When  p  is  a 

non-negative  integer,  i.e.  £>  =  0,1, 2, 3, .  Laguerre  differential  equation  has  as 

solutions  the  associated  Laguerre  polynomials  Llp( x) .  The  first  few  polynomials 
are  listed  in  Table  2. 


Table  2.  Associated  Laguerre  polynomials  From  [9], 

From  equation  (V.25),  the  Laguerre  polynomials  specify  the  constant  £  as 

^-  =  p  so  that  (V.26) 

2  w 
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We  can  now  use  this  result  in  the  equation  (V.20)  in  order  to  find  the  last 
unknown  function  P(z) 
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l  +  l  ilw  °  2p  °  il  i  It  lip 
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We  rewrite  the  last  equation,  using  q(r)  from  (V.22)  and  w2  =  za  +  r/z0 ,  so  that 
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Finally,  integration  from  0  to  r  gives 
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At  this  point,  we  have  found  all  the  unknown  functions  P(r),  q(r)  and  power 
series  Ll  and  we  can  put  them  into  the  trial  solution  (V.7).  After  rearranging 
terms,  we  can  write  down  the  solution  to  the  parabolic  wave  equation  of  the  ( p ,  /) 
mode.  This  solution  is  usually  referred  as  the  ( p ,  /)  Laguerre-Gaussian  mode 
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with  phase 
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(V.  31) 


and  the  Laguerre  polynomials  listed  in  Table  2. 
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C.  LAGUERRE-GAUSSIAN  MODE  ANALYSIS 

1.  Intensity  Plots 

Solving  the  parabolic  wave  equation  in  cylindrical  coordinates,  we  find  a 
family  of  modes,  in  equation  (V.30),  that  are  expressed  in  terms  of  the  product  of 
a  Gaussian  envelope,  a  Laguerre  polynomial  and  a  phase  term.  Each  mode  is  of 
the  (p,l)  order,  where  p  is  the  radial  mode  number  and  /  is  the  angular  mode 
number.  For  this  reason,  the  Laguerre-Gaussian  modes  are  usually  denoted  as 
LG'p  modes. 

Figure  15  shows  the  intensity  profiles  of  various  LG'p  modes  at  r  =  0, 
created  from  simulations  using  C  code.  In  these  plots  there  is  no  phase 
information  displayed,  since  we  are  plotting  the  ap  l  x  a* ,  quantity.  Examining 

these  plots,  we  are  able  to  evaluate  several  features  of  the  Laguerre-Gaussian 
beams. 

We  see  that  the  lowest  order  mode  TG0°  is  just  the  fundamental  Gaussian 
mode,  as  expected.  The  other  modes,  with  p  =  0  and  /  >  0 ,  have  a  simple  ring 
profile,  since  they  have  one  central  node.  As  the  angular  mode  number  / 
increases,  the  central,  near-zero  intensity  region  becomes  larger,  which  is  also 
true  for  modes  with  p  >  0 .  All  modes  with  p  >  0  have  a  profile  in  the  form  of 
concentric  rings.  The  radial  mode  number  p  indicates  the  number  of  nodes  and 
therefore  the  number  of  the  rings.  Thus,  the  total  number  of  dark  nodes  is  p  at 
r*  0  ,  in  addition  to  the  central  node.  These  p  + 1  nodes  come  from  the  zeros  of 
the  Laguerre  polynomials.  We  also  see  that  the  width  of  the  rings  gets  narrower 
for  beams  with  higher  /  and  same  p.  Finally,  as  the  radial  mode  number  p 
increases,  the  effective  width  of  the  beam  gets  bigger. 
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Figure  15.  Laguerre-Gaussian  modes  LG'P  ;  Intensity  plots  (ap[xa’()  in  the 
transverse  plane  x,  y,  at  t  =  0,  where  red  is  most  intense,  while 
blue  is  least  intense. 


In  order  to  see  more  clearly  how  an  LG'p  beam  varies  along  the 

transverse  plane,  we  plot  in  Figure  16  for  two  sample  modes,  the  radial  intensity 
distribution  of  the  Laguerre-Gaussian  beam  across  a  diametric  plane. 
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Figure  16.  Radial  intensity  distribution  for  (a) LG\  and  [b)LG]  beams. 


2.  Plots  of  Real  Part 

Plotting  only  the  intensity  of  a  Laguerre-Gaussian  mode,  we  neglect 
phase.  These  effects  turn  out  to  play  an  important  role  in  the  evolution  of  a 
propagating  beam.  In  Figure  17,  we  plot  the  same  Laguerre-Gaussian  modes  as 
in  Figure  15,  but  this  time  we  include  the  phase  term  of  equation  (V.30)  by 
plotting  the  real  part  squared.  We  could  as  well  plot  the  real  part,  but  in  this  case 
we  would  have  to  plot  negative  values  of  the  optical  field  amplitude  a.  We 
immediately  note  that  in  all  modes,  except  for  the  fundamental  LG°0 ,  radial  nodes 
appear.  The  number  of  nodes  is  proportional  to  the  angular  mode  number  /, 
since  when  r  =  0,  the  phase  (V. 31)  is  just  <j),p(r,&,T)  =  l&. 
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Figure  17.  Laguerre-Gaussian  modes  LG'p  ;  plots  of  the  real  part  squared,  at 

T=  0  . 
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D.  PROPAGATION  OF  LAGUERRE-GAUSSIAN  MODES 

We  have  seen  previously  in  Figure  12  that  a  Gaussian  beam,  which  is  an 
beam,  retains  its  shape  as  it  propagates.  We  find,  by  running  similar 

simulations,  that  this  is  true  for  all  higher-order  modes,  provided  the  laser  beam 
consists  of  only  one  pure  mode. 

1.  Intensity  Plots 

Figure  18  shows  the  results  of  computer  simulations  for  the  evolution  of 
several  Laguerre-Gaussian  modes  as  they  propagate  from  r  =  0  to  r  =  1 .  We 
show  the  intensity  plots  of  the  transverse  cross  section  at  the  beginning  and  at 
the  end  of  the  propagation,  and  the  cross  section  of  the  beam  from  r  =  0  to  r  =  1 . 
It  is  clear  that  as  the  beam  evolves,  its  shape  is  unchanged.  The  central  node  is 
preserved,  and  even  though  the  beam  diverges,  the  multiple  radial  nodes  persist. 
Nevertheless,  we  can  see  clearly  that  the  peak  intensities  of  the  rings  (antinodes) 
decrease  as  a  function  of  the  propagation  distance,  showing  that  the  beam 
spreads  reducing  intensity. 
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Figure  18.  Intensity  plots  of  Laguerre-Gaussian  modes  as  they  propagate. 
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2. 


Plots  of  Real  Part 


In  order  to  complete  our  study  of  the  evolution  of  the  Laguerre-Gaussian 
beams,  we  must  examine  the  phase  structure  of  the  beam  as  it  propagates. 
Simulations  have  shown  that  we  have  to  take  into  account  the  well-known  fact 
that  light  carries  angular  momentum  [10].  For  example,  in  the  case  of  a  circularly 
polarized  beam,  light  caries  an  angular  momentum  that  relates  to  the  spin  of 
individual  photons.  Furthermore,  a  light  beam  can  also  have  orbital  angular 
momentum,  a  feature  that  can  be  deduced  classically  from  Maxwell’s  equations. 
It  has  been  shown  that  for  a  Laguerre-Gaussian  beam,  the  orbital  angular 
momentum  is  well  defined  and  has  a  value  of  Ih  per  photon,  where  /  is  the 
angular  mode  number  and  h  is  the  Planck’s  constant  [1 1], 

The  orbital  angular  momentum  arises  from  the  azimuthal  dependence, 
particularly  from  the  term  el 19  of  the  solution  to  the  wave  equation  (V.30).  We 
will  analyze  this  solution  a  bit  more  by  examining  a  sample  of  Laguerre-Gaussian 
beams.  We  will  plot  the  real  part  of  the  solution  (V.30),  so  as  to  include  the  phase 
and  see  how  it  changes  along  the  propagation  distance.  In  the  following  figures 
we  show  the  surface  plots  of  a  few  Laguerre-Gaussian  modes  at  the  end  of  the 
propagation  distance  (where  r=l).  In  this  way,  we  can  see  how  the  phases 
evolve,  compared  to  their  initial  shapes  at  r  =  0,  presented  in  Figure  17.  For 
each  mode  we  plot  the  real  part  and  the  real  part  squared,  so  that  in  the  latter  we 
avoid  plotting  negative  values,  as  we  did  in  Figure  17. 

Starting  with  the  LG3  mode,  we  put  r  =  1  in  the  solution  (V.30)  and  we 
use  Table  2  to  find  the  associated  Laguerre  polynomial  with  p  =  0  .  Then,  the  real 
part  of  the  solution  can  be  written  explicitly  as 


a0,3  if  ,S,\)  =  a0 


r2z0 

'0 

- 

-r3  e  z°  +1  cos 

r2 

— - - 4  arctan 

+  3,9 

_V+1 

.  (V.32) 
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Figure  19  shows  two  surface  plots  at  time  r  =  1,  for  a  LGl  laser  mode.  In  Figure 
19(a)  we  plot  the  real  part,  as  described  by  equation  (V.32),  while  in  Figure  19(b) 
we  plot  the  real  part  squared.  We  notice  that  a  spiral  pattern  appears  in  the  mode 
structure,  due  to  the  fact  that  the  phase  rotates  as  a  function  of  r  and  the  angle 
&  .  We  have  to  stress  that  the  light  does  not  follow  a  spiral  path;  it  is  just  the 
phase  that  changes  in  a  way  that  it  describes  such  a  spiral  path.  This  path,  in 
three  dimensions  is  somewhat  helical  [10]. 


(a)  (b) 

Figure  19.  Plots  of  (a)  real  part  and  (b)  real  part  squared  of  a  LGl  beam  at  the 
end  of  the  propagation  distance  (  r  =  1 ),  with  ao  =  1  and  zo  =  0.4 . 


Other  examples  of  the  evolution  of  Laguerre-Gaussian  beams,  including 
phase  change,  are  shown  in  Figure  20  and  Figure  21,  where  we  present,  as  in 
Figure  19,  the  surface  plots  of  a  LGl  and  a  LGl  laser  beam  respectively. 


(a) 


(b) 


Figure  20.  Plots  of  (a)  real  part  and  (b)  real  part  squared  of  a  LGl  beam  at  the 
end  of  the  propagation  distance  ( r  =  1 ),  with  ac  =  1  and  z0  =  0.4 . 
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The  solution  to  the  parabolic  wave  equation  for  the  LG\  mode  can  be 
written  in  more  detail  as 


4iz 


a2J(r,(9,1)  =  ao  2 

+1 


2  5 


z  2  +1 


+3/' 


e  z°  +1  cos 


z  2  +1 


-6arctan 


U ; 

,(V.33) 


giving  the  plots  of  Figure  21.  We  observe  again  the  spiral  patterns  due  to  the 
phase  evolution.  We  can  see  more  clearly  how  areal  varies  with  r  in  Figure  22, 

where  we  show  the  radial  distribution  of  areal  along  various  angles. 


(a)  (b) 

Figure  21 .  Plots  of  (a)  real  part  and  (b)  real  part  squared  of  a  LG\  beam  at  the 
end  of  the  propagation  distance  (  r  =  1 ),  with  ao  =  1  and  zo  =  0.4 . 


Analyzing  equation  (V.33),  we  recognize  the  term  e  V+1  as  the  radial 
damping  term.  The  behavior  of  the  oscillating  cosine  term  is  easier  to  understand 
when  we  plot  the  set  of  points  in  r  and  &  where  it  is  maximum  (Figure  23).  We 
see  that  due  to  this  term,  the  spiral  pattern  appears.  These  two  terms,  along  with 
the  term  that  comes  from  the  Laguerre  polynomial,  give  the  beam  profile  of 
Figure  21(a). 
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Figure  22.  Radial  distribution  of  a  LG]  beam  along  (a)  ,9  =  0,  (b)  <9  =  /r/2  and 
(c)  ,9  =  tt/4  . 


Figure  23.  Set  of  points  in  r  and  &  where  the  cosine  term  of  equation  (V.33) 
is  maximum. 
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VI.  HIGH-ORDER  MODES  IN  FREE  ELECTRON  LASER 

SIMULATIONS 


In  previous  chapters,  we  have  studied  the  fundamental  Gaussian  optical 
mode  as  well  as  the  higher-order  modes  in  both  rectangular  and  cylindrical 
coordinates  and  their  propagation  through  a  medium  without  any  sources.  The 
primary  tool  for  our  research  was  the  computer  code,  which  enabled  us  to  solve 
the  parabolic  wave  equation  numerically.  To  analyze  the  operating  conditions  of 
an  FEL  inside  the  undulator,  where  the  optical  beam  coincides  with  the  electron 
beam,  is  more  complex.  There  we  must  solve  both  the  wave  equation  (111.52)  and 
the  pendulum  equation  (111.22)  at  the  same  time.  At  NPS,  the  FEL  research  group 
has  developed  computer  codes  that  solve  these  two  equations  simultaneously, 
and  they  can  be  used  to  simulate  various  working  parameters  for  a  free  electron 
laser  system.  These  are  powerful  tools  to  gain  deeper  understanding  of  FEL 
physics  and  determine  the  importance  of  every  operational  parameter.  Perhaps 
more  importantly,  it  helps  us  to  optimize  FEL  performance. 


A.  ELECTRON  BEAM  TILT 

One  must  first  understand  the  FEL  under  ideal  operating  conditions.  But  in 
order  to  cover  more  realistic  situations,  we  must  study  the  system  performance 
under  non-ideal  conditions.  Vibrations  of  a  mobile  system  for  instance,  might 
cause  the  electron  beam  to  enter  the  undulator  off-axis  or  at  an  angle.  We  use 
the  three-dimensional  simulations  (in  x ,  y  and  r)  developed  in  the  FEL 
research  group  to  investigate  how  sensitive  the  FEL  is  to  off-axis  tilting  and 
shifting.  In  the  simulations,  we  use  dimensionless  transverse  coordinates 
x  =  Xs\kj2L  ,  y  =  y^kjlL  and  the  normalized  time  r  =  cr/L.  Three-dimensional 

simulations  were  used  to  analyze  the  tolerance  of  the  Jefferson  Lab  (Jlab)  FEL  to 
a  tilted  electron  beam.  The  simulations  discussed  below,  are  divided  into  two 
parts. 
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1. 


Gain 


In  the  first  part  of  the  Jlab  simulations,  we  explore  the  gain  response  to  the 
misaligned  electron  beam.  In  other  words,  we  focus  our  attention  on  region  1  of 
Figure  24,  which  is  a  graph  of  a  typical  power  build  up  of  an  FEL.  In  this  region, 
we  have  just  turned  on  the  FEL  so  it  has  not  yet  reached  steady-state.  The  gain 
is  proportional  to  the  slope  of  the  power,  in  this  case. 


Figure  24.  Simulation  output,  showing  the  power  P  versus  number  of  passes 
n  in  an  FEL.  In  the  first  region,  power  is  still  building  up  while  in 
region  2,  FEL  has  reached  its  steady-state. 

We  use  the  basic  input  parameters  shown  in  Table  3,  where  we  simulate 
weak  fields  (ao  =  n)  and  small  number  of  passes  n  so  as  to  be  sure  that  we 
explore  the  power  evolution  in  region  1 .  Also,  the  electron  beam  focus  position  is 
tp  =  0.5  ,  so  the  electron  beam  is  focused  at  the  center  of  the  undulator.  We  thus 

run  our  simulations  by  varying  the  tilt  angle  about  the  center  of  the  undulator. 

We  also  normalize  the  tilt  angle  of  the  electron  beam  to  the  undulator 
length  and  the  optical  wavelength 


3, 


3  3 

y°  ^  y° 

■^k/nL  (426 mrad) 


for  the  Jlab  parameters. 


(VI. 1) 
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Rayleigh  Length 

0.42 

Initial  Optical  Field  Amplitude 

ao 

0.001 

Current  Density  j 

20.0 

Number  of  Undulator  Periods 

N 

30 

Number  of  Passes  n 

30 

Electron  Beam  Focus  Position 

tp 

0.5 

Optical  Waist  Position  tw 

0.5 

Cavity  Quality  Factor  Qn 

17.0 

Table  3.  The  basic  dimensionless  parameters  used  in  the  Jefferson  Lab  FEL 
simulations  for  the  Gain. 

We  vary  the  normalized  tilt  angle  3  from  0  to  4.  To  run  the  simulations,  we 
enter  the  parameters  of  Table  3  and  a  range  for  the  initial  phase  velocity  va .  The 
optimum  value  of  v0  is  determined  in  this  way  and  we  then  look  at  the  output  of 
the  run  with  this  specific  phase  velocity  va .  In  Figure  25,  we  show  a  sample 
output  for  3yo  =  2.0 ,  where  the  optimum  initial  phase  velocity  was  found  to  be 
vo  =  5.5. 

The  dimensionless  input  parameters  are  shown  in  the  upper-right  box  of 
Figure  25.  On  the  upper-left  graph,  we  plot  the  evolution  of  a  slice  of  the  optical 
beam  at  the  outcoupling  mirror,  over  the  number  of  passes  n  ,  along  with  the 
color  scale  for  the  intensity  plots  of  the  optical  field  amplitude.  Next  to  it,  there  is 
a  3-D  representation  of  the  beam  shape  after  the  last  pass,  again  at  the 
outcoupling  mirror.  Below  the  3-D  graph,  we  plot  the  cross-section  of  the  optical 

beam  |a(0,_y,  r)|  (yellow)  and  the  electron  beam  (red)  inside  the  undulator,  at  the 
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last  pass.  Further  down,  there  is  a  graph  of  the  evolution  of  the  phase  velocity 
distribution  over  the  number  of  passes  f(v, n),  and  next  to  it  there  is  a  phase- 
space  plot  of  the  electrons  at  the  end  of  the  undulator.  Below  the  parameter 
table,  there  are  two  plots  which  show  the  development  of  the  optical  power  P(n) 
and  the  gain  G(n) .  In  the  next  row,  there  are  graphs  at  both  the  outcoupling  and 
the  left  mirror  after  the  last  pass,  showing  the  optical  beam  transverse  profile 
(yellow)  compared  to  the  theoretical  Gaussian  (purple),  and  the  optical  phase. 
Also  shown  in  the  same  graphs  is  the  intensity  of  the  light  (using  the  color  scale 
at  the  top)  and  the  extent  of  the  mirrors  (black).  At  the  bottom,  there  is  a  list  of 
the  calculated  results  and  in  the  lower  right  we  plot  the  coefficients  of  the  modes 
|c(m,n)|  that  comprise  the  beam  at  the  outcoupling  mirror  after  the  last  pass.  The 

modal  analysis  is  performed  in  terms  of  the  Hermite-Gaussian  modes.  In  the 
simulation  results  of  Figure  25  we  see  that  the  output  beam  comprises  basically 
from  the  fundamental  mode  (m  =  n  =  0  )  but  simulation  predicts  that  there  will  be 
few  higher-order  modes  in  the  beam  structure  as  well.  This  is  shown  with  the 
plotted  squares  (lighter  blue  color).  Using  this  tool,  we  can  look  at  the  mode 
composition  of  any  arbitrary  wavefront  that  our  simulations  predict  as  the  laser 
output. 

Figure  26  is  a  summary  of  the  results  showing  gain  G  versus  electron 
beam  tilt  angle  9  We  note  that  the  FEL  can  operate  pretty  well  with  an 
imperfect  electron  beam  injection  into  the  undulator,  since  the  gain  falls  to  half  of 
its  peak  value  at  dimensionless  tilt  angle  of  about  9yo  ~  3.5  ,  corresponding  to  an 
actual  tilt  of  about  1.5  rnrad . 
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nx=200,  nt=50,  np=30000,  Wp=8,  seed=7  Thu  Feb  15  11:10:07  2007 

I  I  FEL  WaVEFROHT  EVOLUTION 
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Figure  25.  Output  of  a  3-D  Simulation  for  Jefferson  Lab  FEL. 
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Figure  26.  Summary  of  the  simulation  results  for  the  gain  versus  increasing 
electron  beam  tilt  angle  9 


2.  Extraction 

Having  explored  how  the  FEL  responds  to  a  tilted  electron  beam  during 
the  early  stages  of  its  operation,  we  can  now  explore  the  effects  during  its  steady 
state  phase,  region  2  of  the  graph  in  Figure  24.  So  we  choose  strong  initial 

optical  fields  ao  and  large  number  of  passes  n  ,  as  shown  in  the  Table  4. 

As  a  measure  of  the  performance  of  the  FEL  in  this  case,  we  estimate  the 
extraction  77,  defined  as  the  fractional  energy  transferred  from  the  electron  beam 
to  the  optical  field  during  one  pass  through  the  undulator 

extracted  optical  power 

1  =  .  ~  7  .  7~~T~  - ■  (VI-2) 

initial  electron  beam  power 
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Rayleigh  Length 

0.42 

Initial  Optical  Field  Amplitude 

ao 

10 

Current  Density  j 

20.0 

Number  of  Undulator  Periods 

N 

30 

Number  of  Passes  n 

500 

Electron  Beam  Focus  Position 

0.5 

Optical  Waist  Position  tw 

0.5 

Cavity  Quality  Factor  Qn 

17.0 

Table  4.  The  basic  dimensionless  parameters  used  in  the  Jefferson  Lab  FEL 
simulations  for  the  extraction. 

Running  the  simulations  as  before,  we  first  determine  the  value  of  the 
phase  velocity  v0  for  peak  extraction  77,  and  then  we  look  at  the  output  of  this 
particular  run.  One  sample  output  for  <9V0  =  3.0  is  shown  in  Figure  27,  where  we 
can  verify  from  the  power  evolution  graph  P(n) ,  that  the  FEL  is  in  its  steady-state 
operation  region. 
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nx=200,  nt=50,  np=30000,  T'^p=8 ,  seed=7  Tue  Feb  20  14:23:38  2007 
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Figure  27.  Output  of  a  3-D  Simulation  for  Jefferson  Lab  FEL,  exploring 
extraction. 
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Summarizing  the  results  for  many  tilt  angles  9yo,  we  note  with  some 
surprise  that  the  FEL  exhibits  good  tolerance;  the  extraction  actually  increases 
with  increasing  angles  9yo,  up  to  a  point,  as  shown  in  Figure  28.  This  may  be  due 

to  “overbunching”  of  the  electrons  when  there  is  no  tilt,  reducing  the  extraction  in 
that  case.  The  simulation  predicts  extraction  greater  than  1%  for  normalized 
beam  tilt  9yo  <  3 ,  corresponding  to  an  actual  tilt  about  1.3  mrad . 


Figure  28.  Summary  of  the  simulation  results  for  the  extraction  ?j  versus 
increasing  electron  beam  tilt  angle  9yo . 
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VII.  CONCLUSION 


In  this  thesis,  after  reviewing  the  basic  FEL  and  optical  theory,  we  found 
solutions  to  the  parabolic  wave  equation  that  suggest  the  existence  of  higher- 
order  optical  modes,  working  with  either  rectangular  or  cylindrical  coordinates. 
Motivated  by  the  symmetry  of  the  FEL  system,  we  focused  our  attention  on  the 
solutions  derived  using  cylindrical  coordinates,  which  employed  Laguerre 
polynomials.  Analyzing  the  Laguerre-Gaussian  modes,  we  found  interesting 
results  for  the  beam  shape  and  for  the  beam  evolution  as  well.  The  modal 
analysis  showed  us  that  we  have  to  consider  the  physical  attributes  of  light,  in 
order  to  explain  the  unexpected  behavior  of  the  laser  beam  after  its  propagation. 
Particularly,  we  noticed  spiral  patterns  in  the  beam,  attributed  to  the  orbital 
angular  momentum  of  light.  Finally,  we  demonstrated  the  appearance  of  higher- 
order  modes  in  an  FEL  output,  while  investigating  the  tolerance  of  the  Jefferson 
Lab  FEL  to  electron  beam  tilt. 

Currently,  the  NPS  FEL  research  group  is  analyzing  the  output  wavefronts 
of  computer  simulations  in  terms  of  Hermite-Gaussian  modes.  We  suggest  that  a 
future  thesis  work  would  develop  a  function  that  analyzes  an  arbitrary  wavefront 
from  an  FEL  simulation  in  terms  of  the  Laguerre-Gaussian  modes.  The  use  of 
this  tool  could  make  the  optical  beam  analysis  of  any  simulated  FEL  output  more 
complete,  helping  us  understand  better  the  effects  of  the  investigating  operating 
parameters. 
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